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We consider the non-Markovian Langevin evolution of a dissipative dynamical system in quan- 
tum mechanics in the path integral formalism. After discussing the role of the frequency cutoff for 
the interaction of the system with the heat bath and the kernel and noise correlator that follow 
from the most common choices, we derive an analytic expansion for the exact non-Markovian dis- 
sipation kernel and the corresponding colored noise in the general case that is consistent with the 
fluctuation-dissipation theorem and incorporates systematically non-local corrections. We illustrate 
the modifications to results obtained using the traditional (Markovian) Langevin approach in the 
case of the exponential kernel and analyze the case of the non-Markovian Brownian motion. 

I. INTRODUCTION 

The quest for understanding the evolution of quantum systems under the influence of an environment towards 
equilibrium is ubiquitous in theoretical physics. In particular, the study of quantum-mechanical open systems has 
produced a long list of different techniques and models P, As a matter of fact, most realistic physical systems 
will go through a transient nonequilibrium regime before thermalization is achieved. In the equilibration process, 
interactions with an infinite set of degrees of freedom, which drain energy from the system, usually play a major role. 
This process can in many cases be described stochastically by the use of a Langevin equation [l[ , and was successfully 
applied in a variety of situations. For a comprehensive list of references, see Ref. [2j. 

The microscopic derivation of a Langevin equation in principle yields a dissipation kernel, which encodes memory 
effects or retardation due to the finite reaction time of the environment, and colored noise. These two contributions 
are tied together by the fluctuation-dissipation theorem In a very peculiar limit, namely that of very long times 
compared to the reaction time, one can use the Markovian approximation in which the dissipation kernel is reduced to 
a local term and the noise becomes white. For situations in which such a hierarchy of scales appears naturally, or for 
the cases where one is not interested in the details of transient regimes, this is a very useful simplifying approximation. 

On the other hand, there has been an increasing interest in the nonequilibrium dynamics of phase conversion in a 
number of systems in which transient non-Markovian effects seem to be of relevance. When a precise determination 
of different and yet similar time scales in the process of thermalization must be achieved, one is obliged to include all 
the details in the analysis of the dynamics. In most cases the appropriate framework seems to be that of in- medium 
nonequilibrium field theory J3j, especially in applications to the regimes of nucleation and spinodal decomposition 
after a first-order transition [J] in a myriad of systems, from the formation of the quark-gluon plasma in high-energy 
heavy ion collisions Q to primordial phase transitions in the early universe Q . In the case of ultra-relativistic heavy 
ion collisions, we have shown that memory effects arc important in the determination of the relevant time scales of the 
phase conversion process for the chiral transition [7[ and for the hadronization (confinement) of the deconfined plasma 
as it expands and cools down below the critical temperature 0, and can significantly affect the spinodal explosion 
that presumably occurs [§, [l(| , bringing consequences to the phenomenology @, § • 

Since the structure of memory integrals and colored noise that appears in realistic field-theoretic descriptions of the 
dynamics of phase transitions is often rather complicated 0, HH, El EH, EH EEI , systematic analytic approximations, 
as well as efficient numerical methods, are called for. In this paper, we are interested in the development of analytic 
approximations that can be also useful when coupled with numerical methods. We choose to build our approach in 
the much simpler case of dissipation in quantum mechanics, where all approximations and important scales are under 
control, and where one also finds a wide variety of applications 0- To quote a very relevant, concrete example in 
condensed matter physics, consider the description of localization phenomena in low-dimensional disordered quantum 
systems. The characterization of conductivity properties in a disordered low-dimensional system appears to be well- 
described by a generalized Langevin equation (GLE) with a fully non-Markovian kernel [l||, associated with the 
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predominance of a single frequency in the heat bath [T3|. In this case, one can not disregard a priori non-local effects, 
since they are the essential feature. In this vein, the analysis we present here provides a first step towards a systematic 
semi-analytic means of estimating non-local contributions in a general semiclassical description of quantum dissipative 
systems. Memory kernels also appear in models of financial market data [18l| . In particular, long-range memory seems 
to be present in stochastic processes underlying financial time series which can be originated from market activity 
i.e. the number of trades per unit time [l9j |. In chemical and biological problems Langevin equations with a memory 
kernel are also ubiquitous. For example, the equilibrium fluctuation of the distance between an electron transfer donor 
and acceptor pair within a protein molecule [20| has been shown to undergo subdiffusion and has been modeled by a 
GLE [23. 

The framework that is particularly suited for the integration over degrees of freedom associated with the heat bath, 
and that is most amenable to generalization to field theory is that ofpath integrals I2,[23j]. In particular, it has been 
very successful in applications to the case of the Brownian motion Q HH, H^, [2(1 [27 1 . 

In this paper we study the effects of non-Markovian corrections to the dynamics of a dissipative metastable system 
in quantum mechanics. Starting from the noncquilibrium evolution of a particle coupled linearly to a set of harmonic 
oscillators in the Caldeira-Leggett fashion, we study the effects of the non-local dissipation kernel as well as the colored 
noise that appear in the complete Langevin equation for the particle coordinate in space, q(t). The memory kernel 
has its origin in the Feynman influence functional of the heat bath [28| and is generally too complicated to be treated 
analytically. In the case of field theory, even a numerical analysis is in most cases quite involved. 

To approach the kernel in a simpler, analytic way we develop a systematic expansion in time derivatives of q whose 
convergence is regulated by increasing powers of the frequency cutoff in the distribution of oscillators, Q. Physically, 
above a certain maximum frequency the quantum particle should be "blind" to the bath of oscillators. As we will 
discuss later, one can implement this cutoff in a variety of equivalent ways, in the sense that the only relevant 
parameter to control the time correlation is the width of the distribution. Nevertheless, not all cutoff functional forms 
are allowed if one is concerned with recovering the usual Markovian Langevin dynamics, with a white noise, in the 
limit fit — > 00. Reasonable choices always yield localized kernels, allowing for truncations in the derivative expansion 
that are consistent with the fluctuation-dissipation theorem order by order. 

Since non-Markovian corrections can be very relevant in the description of the dynamics of phase conversion in a 
field-theoretic framework, we try to keep our discussion of the quantum-mechanical case as general as possible, and 
discuss some similarities and differences between the two settings especially regarding the coupling to the environment. 

The paper is organized as follows. Section II presents the quantum-mechanical non-Markovian Langevin equation 
derived within the path-integral framework. We discuss the model adopted for dissipation, based on the Feynman- 
Vernon influence functional, as well as the nonequilibrium dynamics. In Section III we consider some microscopic 
aspects contained in the memory kernel and the noise correlator. In Section IV we introduce and discuss in detail 
the systematic non-local expansion for the memory kernel and the colored noise terms in the generalized Langevin 
equation, describing separately the case of the non-Markovian Brownian motion and dissipative systems subject to 
general potentials. In Section V we illustrate our expansion in the case of an exponential kernel, comparing it to 
well-known exact results for the Markovian Brownian motion and discussing the role of non-Markovian corrections. 
We show that, in this case, the analytic corrections overestimate the non-Markovian effects on the evolution and can 
be used as a rough estimate of their importance before a numerical study. Section VI contains our final remarks. 



II. DISSIPATION AND NON-MARKOVIAN LANGEVIN EQUATION 

We consider a particle, with coordinate q, interacting linearly with a heat bath which introduces dissi pat ion in 
the dynamical evolution. The reservoir is modeled by an infinite set of harmonic oscillators, Rk [H, l29l l3Cl||. We 
assume that oscillators of the bath couple only to the system of interest, in the usual Caldeira-Leggett fashion [3Q], 
and neglect interactions within the reservoir. The hamiltonian is thus given by 

H = H a + H R + Hi (t) , (I) 

where we have defined 

H = ^ + V(q), (2) 



H 



R 



H,(t) = e(t)q^c t R k , (4) 
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and the step function 9(t) simulates the sudden immersion of the system into the reservoir, which is consistent with 
the hypothesis of the nonequilibrium formalism to be used later. 

Using the fact that the interaction is turned on abruptly at t = 0, we assume an initially uncorrelated density 
matrix, with the particle and the bath being in equilibrium independently: 

p(* < 0) = (|0><01), <8» e-^« , (5) 

where T = 1/(3 is the reservoir temperature. 

In order to study the evolution towards equilibrium, we adopt the closed-time path framework [3l[ (for details, see 
Ref. [32]). As is well-known, the method makes use of a complex time contour through which one defines a (complex) 
time ordering. This technique is especially adequate to guarantee the appropriate operator sequence in expressions 
for correlation functions at the expense of duplicating the number of degrees of freedom. In the Keldysh contour, this 
corresponds to a distinction between forward and backward evolution in time, which we denote by indices (+)/(—) 
below. Thus, the generating functional for correlation functions can be computed within the path integral formalism 
in the usual fashion, but with coordinates defined in each of the two branches of the time contour. 

The generating functional is given by 

Z[J + ,J-] = W J [Dq + ][Dq-][DR + k ][DR- k ] 

k 

x expj^y dt(L[q + ,J+,R+}- L[q-,J-,R-])Y 



(6) 



where 

L[q a ,J a ,R a k ] = L[q a ,Rl] + J a q a , (a = ±) . (7) 

Here L [q, Rk] is the lagrangian corresponding to (p}, and J° are auxiliary external currents. 

Since our interest resides exclusively in the evolution of the particle q, we can integrate over the bath variables, 
defining the influence functional [28|: 



*[q+,q-] = ]\([DR + k ][DR7 k ]e^{i( °° dt 

i. J I J-oo 



k 

X 



{L R [q+,R+}+L I [q+ } R+} - L R {q-,R k ]- L:\q-, R^]}^ 



(8) 



The integrals over 

R t e R k 

are quadratic and can be evaluated exactly, yielding: 

■ r + oo 

ln&[q+,q-] = - / dtdt' ]T J a k (t) Gf{t-t') J b k (t'), 

ZJ -°° k;a.b 

(9) 

where a and b are indices for the upper (+) and lower (— ) contours, with G k b (t — t') given by 

Gt + = G>(t-t')9(t-t') + G<(t-t')9(t'-t), 

Gf;~ = G<(t-t'), (10) 
Gk + - G> (t-f), 

G k - = G>(t-t')9(t' - t ) + G<(t-t , )9(t-t'), 
in terms of the retarded and advanced propagators: 

G>(t-t') = — [ e -*-*(*-*') + +e -/Sw* e -*-*(*-*')l (11) 

2mujk 1 — e L J 

G<(t-t<) = -[G>(t-t')]* . (12) 

We are interested in studying the Langevin dynamics in a semiclassical approach. In this vein, we consider the 
evolution of the classical component of q in the presence of thermal and quantum fluctuations. For that purpose, we 
can define a set of more convenient variables by applying a Wigner transform, i.e.: 

(q + ,q-) — (Q,r)=(\(q + +q-),(9 + -q-)) » (13) 
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where Q can be interpreted as a classical solution, while r plays the role of the quantum fluctuations, assumed to be 
small. 

Using Eqs. ([6]) and ([9]), we can express the generating functional as Z = J[DQ][Dr]e lScff , in terms of the following 
effective action: 



S eff = 



dt MfQ-£2(-) 

n-0 





2Jo 



2n+l y( 2 "+ 1 )(Q) 

(2n+l)! 

dt / dt' r(t) ai {t~t')2Q(t') 
Jo 

X_ />OC 

/ dt' r{t)a R {t - t')r{t') . 



where 



a/(t-t') = ]T 



: [ujk{t — t')] coth 



( $Uk_ 
V 2 



2mujk 



sin [wfc(t - f')] • 



In Eq. (fl4]) we also expanded the potential [V(g + ) — V(q )] around Q as follows: 

^ 2 V2/ (2n+l)! ' 

n=0 V ' 

The imaginary exponential of the last term in (|14|) is formally equivalent to 

ypf] exp dt^ dt' t(t)K(t-t')Z(t') 



+ i 



dt r(t)at) , 



(14) 



(15) 



(16) 



(17) 



with X(t — t') = [a^(t — t')] _1 - The variable £(f) enters the equation of motion for Q(t) as a force term. 

The quadratic piece of the exponential above represents the weight for computing averages over £(t), implying 
= an( i (C(^)C(^')) — a R,{t — t'). Assuming that the classical coordinate has a much slower dynamics as 
compared to thermal and quantum fluctuations, we can treat the variable £ stochatically, as a noise term. 

Keeping terms up to 0(r 2 ) and functionally differentiating the action with respect to r, we obtain the equation of 
motion for Q 

(18) 



MQ + V'(Q) -2 dt' ai (t - t') Q(t') = , 
Jo 

which corresponds to the well-known non-local generalization of the traditional Langevin equation. 



III. FREQUENCY CUTOFF, MEMORY KERNEL AND NOISE CORRELATOR 

In order to compute analytically the kernels defined in (|15[) , let us consider the case of a continuous distribution of 
frequencies for the oscillators in the heat bath, p(u>), in the high-temperature limit, T to. 

Physically, there should be a maximum frequency for the interaction of the system with the heat bath. Formally, 
this will play the role of a cutoff scale in our continuous distribution function. One can implement this upper limit in 
different ways. Typical choices are: (i) a sharp cutoff; (ii) an exponential suppression; (hi) a Gaussian distribution; 
(iv) a Lorentzian distribution. The latter provide, as we will see, an exponential kernel, which is the most frequently 
adopted functional form in memory integrals. All these choices have in common the fact of being localized and the 
need of a unique parameter, f2, that characterizes the width of the distribution. In what follows, we consider the 
four cases mentioned above to illustrate possible forms of the memory kernel and which frequency distributions they 
correspond to. For that purpose, we write 

p(uj)c\u) = ^^-f(u;/il), (19) 
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where f(x) is a cutoff function. The particular choice of parameters above is such that all the microscopic information 
concerning the reservoir is encoded in the phenomenological dissipation coefficient r\. 

For each case, we compute the functions an and ai, obtaining results that can be expressed in the general form: 

a R (t-t') = 2r]TA n (t-t'), 

ai{t-t') = ^Ao(t-t'), (20) 

where A$i(i — t') is a representation of the Dirac delta function in the limit SI — > oo. For our choices we find: 

(i) f{u/Q) = 8(1 - w/Q) (Sharp cut): 

»o. 1 sin \Q (t - , \ 

Af l (t-t') = — -. 21) 

° V ' 7T t-f V ' 

(ii) f(u/Cl) = e~ u / n (Exponential cut): 

^g(«-0 = z L JJ* 1 ,,^ )- (22) 



(iii) /(w/O) = e- w / n (Gaussian cut): 



! e -n^(i-t')V4 

^*-*') = ^-^— (23) 



(iv) f(ui/Q) = [{oj/tt) 2 + l] 1 (Lorentzian cut): 



2 

The final form for the equation of motion for the classical coordinate Q is 



AS(t-f) = ?e- n|t -*' 1 - (24) 



MQ + V (Q) + 2 V dt' A n (t - t') Q (*') = ^ (t) , (25) 



o 

where V (Q) is the derivative of the modified particle potential with respect to Q, 

V' (Q) = V (Q) - 2r?A o (0)Q (t) , (26) 

and £ (t) is a colored noise whose two-point correlation function satisfies 

(at)Ht'))='2vTA n (t-t'). (27) 

It is straightforward to show that in the limit Q(t — t') — > oo one obtains the usual Markovian description, with white 
noise. It should be noticed that although the cutoff functions considered above produce slightly different memory 
kernels and noise correlations, the physics is controlled in all these cases solely by the width SI. 

A simple effect that is commonly overlooked in the derivation of a medium-induced Langevin dynamics is the 
modification of the curvature of the original potential brought about by the coupling of the quantum system to the 
heat bath. As can be seen from Eq. (|26p . the correction has always the same functional form, differing only by an 
overall numerical factor of order one characteristic of each cutoff function. 

However, it is customary to ignore this contribution by using the argument that it corresponds to a medium 
renormalization of the bare curvature. Since one is generally interested in comparing results for processes such as 
tunneling or decays in the presence or in the absence of dissipation, one can argue that one should keep the original 
potential fixed 0, [3(| ■ This can be attained by introducing appropriate counterterms in the lagrangian. 

The procedure of treating differently the modifications on the original potential and the effect of dissipation on 
the movement of a particle in the framework of quantum mechanics is usually well justified [2], |30(. In this case, the 
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degrees of freedom of the reservoir are totally independent of those of the system of interest by construction, so that 
one can use the strategy above to pinpoint specific effects of dissipation on, for instance, a decay rate 1 . 

In the case of a field theory, though, dissipation is commonly introduced by isolating what can be seen as a classical 
part, or a condensate, of the field from its other modes. The latter play the role of heat bath modes. Nevertheless, 
they can not be considered as an external medium in which one will embed the system of interest, now represented 
by a specific mode of the field. System and reservoir are necessarily entangled. Therefore, consistency implies that 
every medium correction should be taken into account. For our choice of frequency distribution, the correction to the 
potential has the form AV = — ?yAo(0)<5 2 . Since the correction AV is proportional to fl, it is clear that it contains 
an ultraviolet divergence, and also that for large enough values of f2 the sign of the curvature may be inverted. The 
limit of very large f2, in spite of yielding the Markovian regime, is physically meaningless regarding the interaction 
of the system of interest with the oscillators of the bath. In reality, this interaction will be significant only within a 
finite window in the spectrum of frequencies of the oscillators. 

With the results above, one has all the ingredients to study the semiclassical evolution of a particle subject to a 
dissipative medium, in the presence of an arbitrary potential. Given the complexity of the memory kernel, one is 
usually forced to make use of numerical techniques in order to obtain exact results. However, as will be shown in the 
next section, one can also resort to a convenient derivative expansion, which introduces non-Markovian effects order 
by order. In this way, one can proceed further with analytic steps and eventually simplify the numerical simulations. 
As discussed previously, this could be especially useful in the context of quantum field theory. 



IV. SYSTEMATIC NON-LOCAL EXPANSION 



A memory integral such as the one in (|25p is, in principle, always present in a realistic Langevin equation. Due to its 
highly non-local nature, it usually has a very complicated structure to be treated analytically, and even numerically, 
especially if a field-theory approach is needed. However, in the simpler case of quantum mechanics, one is able to 
derive a convenient series expansion in powers of I/O, which includes higher non-Markovian effects systematically. As 
will be shown below, the first term of this series yields the usual Markovian approximation to the Langevin evolution 
when the appropriate limit is taken. 

The method resorts to the assumption of a hierarchy of relevant time scales. In addition to the separation which 
occurs naturally in the high-temperature limit, as discussed before, one has to consider the limit of large times, so 
that we end up with 

t > t coU > i , (28) 

where t co u = I/O, is the collision time. 

The limit in which fit is strictly infinite yields a local dissipation term ~ Q(t). The non-Markovian contributions 
arise from the assumption that fit is still large, but finite. As one goes down in fit one should, in principle, include 
more and more non-local corrections in the series for the memory integral. 



A. Non-Markovian Brownian motion 



Before the systematic incorporation of non-Markovian corrections in a general dissipative system, let us consider 
the description of Brownian motion, defined by V (Q) = 0. In this case the method is most easily implemented in 
the Laplace space. It is convenient to introduce a dimensionless (small) expansion parameter 5 = rj/Mfl and work 
with a dimensionless time variable r = (r]/M)t. In terms of these, Eq. ([23]) for the case of Brownian motion can be 
written as 



d 2 Q(r) 



dr 2 

where £ = (M/ij 2 )^ and A$(t — r') is related to An(i — t') as 



2^ T rfr'A 5 (r-r')^P = C(r), (29) 



A s (t-t') = — Ao(t-t'), (30) 
V 



1 A remarkable exception is found in Ref. [3^1 . where the author applies instanton techniques to deal with phonon modes in path integrals. 
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so that the noise correlation function is given by 



MT 

(C(r)C(r'))=2 — A s (t-t'). 



(31) 



Deno ting by s the Laplace-conjugated variable to r and by f(s) the Laplace transform of a time-dependent function 
f(r) M 

(32) 



(33) 



f(a) = C[f(r)] = / dr e~" f{r), 
Jo 

one can solve (|2"91 in the Laplace space algebraically 

Q{s)=g(s) v + Q Q s + 2Q A s (s) + ((s) 
where Q Q = Q(t = 0), v = Q(t = 0), and 

g(s) = 



s + 2A 5 (s) 



(34) 



Up to here, all results are exact and for some kernels the inversion of the Laplace-transformed equation can be 
calculated in closed form without approximations - as is the case for the exponential kernel (|24p . Now, for a kernel 
of arbitrary form, As (s) can be written generically as 



A s (s) = ±[l + R(5s)}, 



(35) 



where i?(0) = 0, since it must reproduce the Markovian limit for S — > 0. Therefore, one can write for g(s) the following 
power series expansion 



i 



s(s + l) s(s + l)^-f 

v J v ' n— 1 



oo / 1 x n oo dI'^ 



(36) 



fc = l 



where the coefficients RqI are the fc-th derivatives of [R(Ss)] n evaluated at 5s — and contain all the dependence on 
the specific form of the kernel As(t — t')\ 



~g{s) 



s + 2A s (s) 



s(s 



^i) { 1+ li ( ~ ir (jtt)" [r({8)1 "} 



(37) 



The power series expansion of [R(5s)) n is simply 



R 



k=l 



(38) 



Replacing the expansion (|36| in (|29p . one obtains an expression for Q(s) that can be easily Laplace-inverted to obtain 
Q(t). To simplify the presentation we consider Q = and vq = 0, so that 



Q(s) 



oo oo , 

7 v ' fc=l |_«=1 



s(a + l) w s(s + l)^! 



(6s) 



(39) 



The inversion gives the general solution for the semiclassical (non-Markovian) evolution of the particle in the harmonic 
heat bath: 



Q(r) = / dr'C(r') 



1 - e 



-(r-r') 



°° (_ i\n fT 

Y, 6 -ir<l dr'C(r')F^(r-r% 

k,n—l 



(40) 
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with 



F (n,fc)( T ) = c -l 



(s + 1)" +1 



1 d^- 1 ) 



i\ dr( k 



(41) 



This is our main result for the general non-Markovian Brownian motion. The first term in Eq. (|40p corresponds to 
the Markovian approximation while higher orders in 5 = Mjrfil are associated with corrections due to non-locality. 

(n) 

Notice that this expression is valid for any memory kernel with well-defined derivatives R - such as the kernels 
in (|2"Tj) - (f2~4")) . We also note that the same method presented above can be applied straightforwardly to the case of a 
harmonic potential, V(Q) ~ Q 2 , since the associated Langevin equation (|25|) describing this system is still linear. For 
a generic anharmonic potential the Laplace transformation method used above is not useful and a different strategy 
must be employed. 



B. General external potential V(Q) 

In what follows, we consider a generic external potential V(Q) and assume that our system will eventually thermalizc, 
so that the particle will end up at a stable minimum of the potential. In this case, we can assume that Q(t') is bounded 
within the interval [0, t], and Q(t — > oo) — > 0. Given this condition, and the fact that 

lim A n (*-*')-»<*(*-*') , (42) 
n — >oo 

the memory integral in (|25[) has support in the vicinity of t' = t for £1 very large but still finite. Defining the 
dimensionless variable x = Q (t — t'), one can expand the memory kernel in p5|) in a Taylor series around x — and 
obtain 

f dt' A n (t - t') Q (f ) = £ ^Q(" +1 ) (t) , (43) 
where we have defined the integral coefficients 

Jniy) = —J dxx n ' (44) 

One should notice that, since the only dependence on fl of the kernel written as above is An(x/f2) ~ ft, the combination 
Afi(x/Q,)/tt is independent of f2, and we have a well-defined series in inverse powers of ft for the memory integral 2 . 

To obtain the corresponding expansion for the two-point correlation function of the colored noise that satisfies the 
fluctuation-dissipation theorem order by order, we can recast the expansion above as an expansion for the kernel by 
using identities involving the derivatives of Q and derivatives of the delta function. Namely, we can write 



f dt' An(t-0 Q(t') = 2 f dt' 
Jo Jo 



^ (-l)"J n (Ot) cT 

.n=0 



Q(t) , (45) 



so that the expansion for the noise correlation function reads 



(m(t')) - 4r,T£ LJ)y^l£-S(t - f ) , (46) 

n—0 

an expression that clearly shows the increase in non-locality as higher-order terms are needed for a given memory 
kernel. The presence of the derivatives, as also happens in some typical stochastic evolution equations such as the 
Cahn-Hilliard equation for conserved order parameters [1, Hi| , naturally suggests that the problem should be solved 
in the Fourier space. 



2 Although f2 is still present in the upper limit of the integral in J n , the integral is dominated by x ss given the localized behavior of 
the kernel Aq. 
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The equation of motion (|25p is, then, equivalent to 

MQ + V' (Q) + 2 V £ ^^Q<" +1 > (t) = f (*) , (47) 

which reduces to the traditional Langevin equation with white noise in the limit fit — > oo, since J n (oo) — > S n o/2 and the 
correlator in ([27]) tends to 2?7T<5 (t — t'), consistently with the fluctuation-dissipation theorem. When one incorporates 
non-local corrections from the expansion of the memory kernel to a given order, one should also incorporate corrections 
to the same order in the noise correlator. In this way, the fluctuation-dissipation theorem is satisfied order by order. 

Inspection of (|47[) shows that terms containing higher-order non-local corrections in time derivatives of Q, and 
correspondingly in time derivatives of the delta function for the noise correlator, are strongly suppressed by increasing 
powers of I/O. This allows for cutting the series at a given value of n. Including just the first two non-local 
corrections 3 , i.e. going up to n = 2, we obtain: 

M(t) Q + V'(Q)+r ll (t) Q + va(t)Q=Z(t), (48) 

where V' (Q) satisfies (f26f and £ (t) satisfies (|46|) with the sum cut at n = 2, and we have defined 

M(t) = M+^J^flt), (49) 
77! (t) = 2 V J Q (Qt), (50) 
%(*) = §^(Oi). (51) 

One should notice that the first non-local corrections not only yield new terms but also modify the ones which were 
already present in the Markovian limit. In general, the mass and the dissipation coefficient acquire time-dependent 
corrections that become constant after a transient period t > 1/Vt (see below). 

V. APPLICATION TO THE CASE OF AN EXPONENTIAL KERNEL 

To illustrate the method, we consider the case of an exponential kernel, which corresponds to a Lorentzian distribu- 
tion of frequencies as discussed previously. The choice of this kernel is motivated by its frequent use in the literature 
and by the fact that it can be solved exactly. It corresponds to the so-called Ornstein-Uhlenbeck process [l|, [24| , and 
one can solve it exactly by converting the non-local problem into a set of Markovian equations. These processes can 
also be used to test the reliability of numerical codes to solve non-Markovian stochastic evolution equations [36|, l37[ . 

In this case, one can easily compute the first coefficients J„, obtaining: 

Mv) = \ [1 - e~ V ] , (52) 
•h{y) = -\[l-e-y (1 + y)] , (53) 



My) -I 



1 



1/(1/ + 2) 



(54) 



2 

It is clear that, except for exponentially suppressed corrections, the coefficients will be all constant, and given by 

Mv) = My) = -Ji{y) = \, (55) 



3 In the cases of sharp and exponential cutoff functions, the integral defining J n (y) actually diverges after some value of n for large y. As 
stated before, we assume that the system eventually thermalizes, so that the corresponding time derivatives of Q vanish fast enough, 
taming each term of the series expansion dynamically. In the more physically sensible (smooth) frequency cuts, such as the Gaussian 
and Lorentzian cases, however, the functions J n are always finite. 
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so that the approximate non-local Langevin equation, to order n = 2, is approximately given by 

V 



M — 



Q + v' (Q) + vQ + -^ Q=Ht) 



(56) 



an equation that is, of course, only valid in the limit of fl very large as discussed before. The noise two-point correlation 
function to this order is given by 



mm) = w 



5(t-t') + ±5>(t-t') + ±6"(t-t') 



(57) 



Since it is an initial value problem, we can solve this equation using the method of Laplace transforms [34j . For 
simplicity, we consider the case of Brownian motion 4 (V (Q) = 0) so that we deal with a linear equation, and the 
following initial conditions: Q(t = 0) = 0, Q(t = 0) = Vq, Q(t = 0) = 0. It is convenient to work with the same 
dimensionless time variable as before, r = (rj/M)t = t/a, and its Laplace conjugate s, as well as with the previous 
(small) dimensionless quantity S = r)/MQ. Converting the Langevin equation to the Laplace space, we can solve it 
algebraically obtaining 



Q(s) = g(s) 



(1 - S)v + S 2 v Q s + -i(s) 



V 



where tilde variables are Laplace transforms of the original ones and we have defined the function 



s [1 + (1 - 8)s + S 2 s 2 } ' 
We can now use the convolution theorem for Laplace transforms to write Q(t) as 



Q(r) = - f dr'ar')g(T-r'), 
V Jo 



with g(r) given by 



g (r) = l + _L [(A - 1 + 2S)e-^ 1+A ^ 2S + (A + 1 - 25)e-^- A) ' 2S 



(58) 



(59) 



(60) 



(61) 



where A = (1 — AS) 1 / 2 . Since S is assumed to be small from the outset and kept to second order in our expansion, 
it is simpler to first expand Q(s) in powers of 6, then perform the inverse Laplace transforms. We show the explicit 
result for vq — 0, for compactness 5 : 



Q(r) 



a 



dr 



V Jo 



'£(/) { (l - e~ (T ~ T,) ) + <5e-( T - T ') (r - t') 
(r-r'r 



2(t - t') - 1 



0{8*). 



(62) 



Using the fact that (£(t)) = 0, it follows that the average classical position vanishes, (Q(t)) — 0. Quantities such 
as (Q 2 (t)} can be expressed as 



(Q 2 {r)) = (Q{r)) 2 + / dr'dr"g{r - r')g{r t")<£(t'X(t")> , 



(63) 



4 For a non- vanishing potential, the procedure would be the same discussed here, but involving more complicated Laplace transforms and 
convolutions. In any case, always integrals to be computed once instead of in each time step of evolution. This method is, of course, 
reasonable only if one does not need to include many higher derivatives to account for the deviation from the Markovian regime. 

5 In textbook treatments of the Brownian motion, one usually assumes a Gaussian distribution of initial velocities [H l24ll . In the limit of 
early times, its mean-square average is responsible for the well-known behavior (Q 2 ) ~ (^q)* 2 . Since we fixed vq = for simplicity and 
compactness of analytic expressions used basically to illustrate the non-Markovian expansion we propose, this limit will not appear in 
our final result for (Q 2 ). Nevertheless, it can be straightforwardly incorporated to reproduce the results of Ref. in the Markovian 
limit. 
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and we know the noise correlator as an expansion in derivatives of the delta function 5(t — t'), Eq. (157)) . Since we 
have already an expansion in S of Q(r), we can compute (Q 2 (t)} directly, obtaining: 



(Q 2 )(r) 



2MT 



V 

,-2t 



[2 - 5(3 
1 

2 + 



2t) 



5 2 (1 



3r 



5(1 



- — (l + 10r - 2r 2 

4 v 



0(<5 3 ). 



In the limit of t > 1, i.e. £ 3> M/tj, we obtain the usual late-time diffusion behavior 



(Q 2 )lateW 



2T 



(64) 



(65) 



whereas for 5 -C r <C 1, i.e. l/fi <C t *C M/ry, we obtain 



T 



(g 2 )carlyW«^t 2 . 



(66) 



ignoring corrections ~ 0(5 3 , t 3 ). As expected, we have essentially two regimes separated in time by the scale a = M/rj, 
which measures the relative importance between the second and first derivatives in the Langevin equation as well as 
the time required to erase the memory of the initial velocity. 

Had we kept a non-zero initial velocity, we would find the well-known early-time behavior Q 2 ss v 2 t 2 also in the 
non-Markovian limit. Therefore, it is clear from Eq. (f6"6"| that the effect of non-Markovian corrections at early times 
is to modify the initial velocity (or average velocity, if one assumes, as customary, a Gaussian distribution of initial 
velocities) of the Brownian particle. Recall that, due to the equipartition theorem, T/M represents the average of the 
square of the velocity of the particle. 

In Fig.[T]we show results for the normalized average quadratic dispersion of the quantum particle (Q 2 )(t) / (2MT frf) 
in the Markovian and non-Markovian cases, with 5 = 0.5 in the latter. The comparison between the exact non- 
Markovian solution and the non-Markovian results up to 0(5), 0(5 2 ), and 0(5 3 ) 6 clearly illustrates the convergence 
of the expansion. 



H 



A 

N 

o 

v 1 



0( 5 ) ; 6 = 0.5 

c-2 



0( 8 ) ; 8 

r.3 



= 0.5 

0( 5 J ) ; 5 = 0.5 
Exact ; 8 = 0.5 
Markov (8 = 0) 




FIG. 1: Normalized average quadratic dispersion of the quantum particle (Q 2 ) as a function of the dimensionless time r = 
(rj/M)t for 5 — 0.5. The usual Markovian result (black, dotted line) is compared to the behavior in the presence of non-local 
corrections up to order ~ 5, ~ 5 2 and ~ <5 3 (dashed lines, top to bottom, respectively), as well as with the exact non-Markovian 
Brownian evolution (magenta, solid line). 



The calculation of the full 0(<5 3 ) non-Markovian result includes the n = 3 term, T]4,(t) Q, in our systematic expansion, being a 
straighforward computation. 
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Notice also that the exact result lies between the Markovian and the approximate non-Mar kovian results. Such 
a counterintuitive behavior can be interpreted as a consequence of the effectively diminished mass that guides the 
truncated non-Markovian evolution, especially at early times: M = M + 2r]Ji/fl. The effective inertia of the non- 
Markovian Brownian particle described by the lowest non-trivial order of our truncated systematic expansion is clearly 
smaller than in the corresponding Markovian approximation, yielding therefore larger quadratic dispersions. Higher- 
order non-Markovian terms bring about higher derivatives in the evolution equation that contribute eventually 
to the quantitative correction of this trend, approaching the exact result. It is interesting to note that this is a general 
feature of the 0(d) result, regardless of which explicit memory kernel is being considered. Since the integral J\ is 
always negative (cf. Eq. (|44]) ). the effective mass M = M + 2r)Ji/tt is always smaller than the original one. Provided 
that the 0(d) non- white noise does not compete with the reduced effective inertia, the average quadratic dispersion 
will be increased in the truncated non-Markovian case. For the case of an exponential kernel, Fig. [1] shows that the 
first non-trivial order of our truncated systematic expansion is in fact an overestimate of the Non-Markovian effects 
in the quadratic dispersion. 

Fig. [2] displays the late-time behavior of (Q 2 ) for different values of S: (local), 0.1 and 0.3 (non-local). It illustrates 
the modifications brought about by the first non-trivial corrections in our expansion for reasonable (not too large) 
values of <5. The normalization is given by the Markovian late-time limit, Eq. (|65p . All curves flatten out for large 
enough times, as expected, with the late-time behavior being exactly the same. 




FIG. 2: (Q 2 ) normalized by the Markovian late-time behavior, Eq. (|65[1 . as a function of the dimensionless time r = (rj/M)t. 
The Markovian approximation (S = 0) is compared to the exact and truncated (up to 0(5 3 )) non-Markovian results for different 
values of 8 — rj/MSl. 



Using the systematic non-local expansion presented in the previous section, one can then estimate analytically the 
deviation from a Markovian behavior of the early-time dynamics of a given system straightforwardly. For instance, 
this can be used as a prescription for evaluating the importance of transient effects generated by non-Markovian 
corrections on the dynamics of the system under investigation to decide whether it is worth the effort of incorporating 
them in a full numerical simulation. Conversely, it can be used to incorporate memory effects order by order in 
5 = rj /Mil without the need of performing a full memory integral in each time step of a simulation. 



VI. CONCLUSIONS AND OUTLOOK 



The complex structures of memory kernels, which usually can be cast into non-trivial dissipation and noise terms 
in a Langevin evolution equation, can not be disregarded even as a first approximation in several systems of interest 
in condensed matter, as well as in nuclear and particle physics under extreme conditions. Subdiffusion problems in 
chemical and biological systems and non-Abelian plasmas are just a few examples of this wide class of phenomena. In 
most cases, however, the incorporation of non-local effects is either overlooked or performed uniquely via numerical 
computation. The latter usually demands heavy computer power and has the caveat of not providing the type of 
insight that one can extract from analytic calculations, even within a simple framework. 
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In this paper we presented a systematic semi-analytic method that can be combined effectively with numerical 
techniques to tackle the non-Mar kovian Langevin evolution of a dissipative dynamical system in quantum mechanics. 
Making use of the Schwinger-Keldysh formalism within the path integral framework, we discussed the physics content 
of the frequency cutoff for the interaction of the system with the heat bath, which is commonly taken for granted in 
the usual Caldeira-Leggett description of open systems in quantum mechanics. We considered the kernel and noise 
correlator that follow from the most common choices, and also some less frequent, and derived an analytic expansion 
for the exact non-Markovian dissipation kernel and the corresponding colored noise for a general potential that is 
consistent with the fluctuation-dissipation theorem and incorporates systematically non-local corrections. 

Using the non-Markovian Brownian motion with an exponential kernel, we illustrated how to implement the method 
and showed its consistency in recovering the relevant and well-known limits. This technique can be applied to any 
kernel that is amenable to a series expansion of the sort performed in this paper, which represents a very wide class 
of possibilities for applications. We believe that a combination of the semi- analytic treatment proposed in this paper 
with regular numerical methods can be of use in the study of several physical systems, at least providing shortcuts 
for numerical calculations and reasonable estimates for the relevance of non-local corrections given the typical scales 
of the problem under consideration. 

Although the Laplace method is much more limited in the case of nonlinear equations, as will be the case once one 
goes beyond the harmonic approximation for the external potential, it is of great power for the Brownian problem 
with memory. For an external potential of generic form one can avoid the Laplace method at the cost of introducing 
higher order time derivatives. With the method presented in this work, one can address pragmatically several physical 
systems that exhibit this behavior, as was discussed in the introduction. The extension to quantum field theory is 
certainly not trivial, but possible. 
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